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Abstract: We compute the matching relation for the strong coupling constant within the 
framework of QCD up to four-loop order. This allows a consistent five-loop running (once the 
(3 function is available to this order) taking into account threshold effects. As a side product we 
obtain the effective coupling of a Higgs boson to gluons with five-loop accuracy. 
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1. Introduction 

The strong coupling constant, Og, constitutes a fundamental parameter in the Standard Model 
and thus its precise numerical value is very important for many physical predictions. An inter¬ 
esting property of Ug is its scale dependence, in particular its strong rise for low and its small 
value for high energies which make perturbative calculations within the framework of QCD pos¬ 
sible. The scale dependence is governed by the (3 function. However, in order to relate ag at 
two different scales it is also necessary to incorporate threshold effects of heavy quarks which is 
achieved with the help of the so-called matching or decoupling relations. Thus, when specifying 
ag it is necessary to indicate next to the scale also the number of active flavours. In this paper 
we evaluate the decoupling relations to four-loop accuracy. This makes it possible to perform a 
consistent running of the strong coupling evaluated at a low scale, like, e.g., the mass of the r 
lepton, to a high scale like the Z boson mass — once the five-loop (3 function is available. 

Many different techniques have been developed and applied to various classes of Feynman 
diagrams. The complexity increases both with the number of legs and the number of loops. As 
far as the application of multi-loop diagrams to physical processes is concerned the current limit 
are four-loop single-scale Feynman diagrams, where either all internal particles are massless and 
one external momentum flows through the diagram (see, e.g.. Ref. [1] for a recent publication), or 
all external momenta are zero and besides massless lines there are also particles with a common 
mass M. The latter case has been developed in Refs. [2,3] and first applications can be found 
in Refs. [4,5]. In this paper we consider a further very important application: the four-loop 
contribution to the matching or decoupling relation for the strong coupling. 

The paper is organized as follows: In the next Section we define the decoupling constants 
and the theoretical framework of our calculation. In Section ^ we present analytical results and 
discuss the numerical consequences. In Section ^ the connection of the decoupling constant to 




the coupling of a Higgs boson to two gluons is explained and the corresponding coupling strength 
is evaluated to five-loop order. Finally, we conclude in Section In the Appendix we present 
the result for the decoupling constant parameterized in terms of the on-shell heavy quark mass. 


2. Theoretical framework 


We consider QCD with Uf active quark flavours. Furthermore it is assumed that ni quarks are 
massless and Uh quarks are massive, i.e. we have nj = ni + rih. In practice one often has Uh = 1, 
however, it is convenient to keep a generic label for the massive quarks. 

The decoupling relations relate quantities in the full and effective theory where the latter is 
defined through the Lagrangian C! given by 


£' (ff°, m°, c°’“; C°) = ( 5 °', 


( 2 . 1 ) 


tpq, and c“ are the fermion, gluon and ghost fields, respectively, niq are the quark masses, ^ 
is the gauge parameter, and as = < 7 ^/(dvr) is the strong coupling constant. is the usual 

QCD Lagrange density and the effective n;-flavour quantities are marked by a prime. Eq. o 
states that the Lagrangian in the effective theory has the same form as the original one with 
rescaled helds, masses and coupling. It is convenient to define the decoupling constants Q in the 
bare theory through 
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In a next step the renormalized quantities are obtained by the usual renormalization pro¬ 
cedure introduced by the multiplicative renormalization constants through [ 6 ] 


9^s=9‘^Zg9s, m'^ = Zmmq, - 1 = ^ 3(1 - 1 ), 

= = c°’“ = V^c“. (2.3) 


Combining Eqs. (|2.2|) and (|2.3|) leads to renormalized decoupling constants, e.g. 


C9 


^9 ^0 


Cs 



Cs 



(2.4) 


Note that since we are interested in the four-loop results for Q the corresponding renormalization 
constants have to be known with the same accuracy. In Ref. [7] the results up to four-loop order 
have nicely been summarized (see also Refs. [8,9]). 

Due to the well-known Ward identities [ 6 ] there are several ways to compute the renormal¬ 
ization constant for the strong coupling, Zg. A convenient relation, which has the advantage that 
due to the appearance of renormalization constants involving ghosts less diagrams contribute, 
is given by 
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where Z\ is the renormalization constant of the ghost-gluon vertex ggGcc. The same is true for 
the corresponding equation for the decoupling constant, such that one can use the relation 
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where denotes the decoupling constant for the ghost-gluon vertex. Alternatively, one can use 
the renormalized objects ^ 3 ; Cs from Eq. (2.4) as well as Ci = §^Ci ^.nd then obtain Q from the 


renormalized version of Eq. 

In Refs. [10,11] formulae for the bare decoupling constants are derived which relate the 
n-loop decoupling constants to n-loop vacuum integrals. In particular, one has 

C3° = i + n^'^(o), 

C 3 ° = i + nf(o), 

C? = i + r^i(o,o), (2.7) 

where nc(p^) and nc(p^) are the gluon and ghost vacuum polarizations, respectively, and the 
superscript h denotes the so-called hard part which survives after setting the external momentum 
to zero. Specifically, nc(p^) and Tldp^) are related to the gluon and ghost propagators through 
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respectively, while ^gcc(p^^) defined through the one-particle-irreducible (IPI) part of the 
amputated Gcc Green function as 

iPi 


(2.9) 


^GcciP^ A:)] -|- other colour structures| , 


where p and k are the outgoing four-momenta of c and G, respectively, and are the structure 
constants of the QCD gauge group. Sample four-loop diagrams for each line of Eq. (2.7) are 
shown in Eig. ||(a)-(c). 

From Eqs. dH), <M) and (^)^) it becomes clear that for the calculation of (g four-loop 
vacuum diagrams are needed. Currently the only practical method to express an arbitrary 
four-loop vacuum integral in terms of a small set of master integrals is based on the algorithm 
developed in Ref. [12]. The application to four-loop bubbles has been discussed in Ref. [2]. 
First physical results deal with moments of the photon polarization function [4] and the singlet 
contribution to the electroweak p parameter [5]. The essence of the Laporta algorithm [12] is 
the generation of large tables containing relations between arbitrary integrals and the so-called 
master integrals. For the calculation at hand the tables have a size of about 8 GB and contain 
6 million equations. 

The master integrals needed for the evaluation of (g have been computed in Ref. [13], where, 
however, some of the higher order coefficients in e could only be determined numerically. 
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Figure 1: Sample diagrams for the gluon (a) and ghost (b) propagator and the ghost-gluon vertex (c). 
In (d) the lowest-order diagram is shown mediating the Higgs-gluon coupling in the Standard Model and 


(e) shows an example for a five-loop diagram contributing to the result in Eq. (4.4). 


3. Running and deconpling for Og 


Whereas at three-loop level of the order of 1000 diagrams have to be considered, at four loops 
there are almost 20000 diagrams which contribute to the gluon and ghost propagators and 
the ghost-gluon vertex. They are generated with the program QGRAF [14]. With the help of the 
packages q2e and exp [15,16] the topologies and notation are adopted to the program performing 
the reduction of the four-loop vacuum diagrams [2]. As an output we obtain the bare four-loop 
results as a linear combination of several master integrals. All of them have been computed in 
Ref. [13]. 

Since at four-loop order the renormalization is quite non-trivial, let us in the following briefly 
describe the procedure necessary to arrive at a finite result. It is convenient to build in a first step 
the sum of the bare contributions to C3 , C3 3 -nd and combine them immediately to (g according 
to Eq. ( |2.6|) . Already at this point the gauge parameter, which for the individual pieces starts 
to appear at three-loop order, drops out and hence spares us from renormalizing Let us 
mention that due to the complexity of the intermediate expressions, the four-loop diagrams 
have been evaluated for Feynman gauge, whereas the lower-order diagrams were computed for 
general 

In a next step it is convenient to renormalize the parameters a* = gg/(47r) and rrifi applying 
the usual multiplicative renormalization (cf. Eq. (|2.3| )). The corresponding counterterms have 
to be known up to the three-loop order. At this point one has to apply Eq. ( p.4| ) which requires 
the ratio ZgjZ'g up to four-loop order. In order to evaluate this ratio one has to remember that 
Z'g is defined in the effective theory and thus depends on and n; whereas Zg depends on Og 
and (n/ -|- Uh)- Thus it is necessary to use (g up to three-loop level in order to transform a'g to 
ag where due to the presence of the divergences in Zg also higher-order terms in e of Cg have to 
be taken into account. 

Finally one arrives at the following finite result for {Cg)‘^ which for Nc = 3 and = 1 is 
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(3.1) 


where the heavy quark mass m^ is renormalized in the MS scheme at the scale jj.. The corre¬ 
sponding expression for the on-shell mass is given in Appendix]^ In Eq. (3.1), C(n) is Riemann’s 
zeta function and a„ = Li„(l/2) = 1/(2^A:”). The constant Xq, which is the leading coef¬ 

ficient of a certain finite four-loop master integral, is only known numerically with the value [13] 


Xo = -M.808879546208334741426364595086952090. 


(3.2) 


Interestingly, in principle the number of numerical coefficients occurring in Eq. (P) should be 
three. One relation among them can be established through the separate renormalization of the 
ghost propagator while a further constant has become available recently in analytical form [17]. 
Thus one remains with one coefficient which is only known numerically. 

Inserting numerical values into Eq. (PD one obtains 
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It is interesting to note that the n;-independent four-loop coefficient is relatively big as com¬ 
pared to the corresponding constants at lower loop-order. However, for the interesting values 
ni = (3,4, 5) one observes a big cancellation leading to a well-defined perturbative series with 
coefficients (-0.4288,-1-0.7790,-1-1.9428) in front of (a^/vr)'^. 

We are now in a position to study the numerical impact of our result. As an example we 
consider the evaluation of a^s\Mz) from i.e. we apply our formalism to the crossing 

of the bottom quark threshold with ni = 4. In general one assumes that the value of the scale 
fj,b, where the matching has to be performed, is of order mb- However, it is not determined 
by theory. Thus this uncertainty contributes significantly to the error of physical predictions. 
On general grounds one expects that while including higher order perturbative corrections the 
relation between {Mr) and ai^\Mz) becomes insensitive to the choice of the matching scale. 
This has been demonstrated in Refs. [10,18] for the three- and four-loop evolution, respectively. 
In the following we want to extend the analysis to five loops. 

The procedure is as follows. In a first step we calculate (/xf,) by exactly integrating the 
equation 


dfi'^ vr 
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with the initial condition {Mr) = 0.36. Afterwards a^J^\fjLb) is obtained from the renormal¬ 
ized version of the first equation in (^) where we use Cg parameterized in terms of the on-shell 
mass (cf. Eq. (A_T)) Mb = 4.7 GeV. Finally, we compute a^p{Mz) using again Eq. (^^). For 
consistency, z-loop evolution must be accompanied by {i — l)-loop matching, i.e. if we omit 
terms of 0{al'^'^) on the right-hand side of Eq. (|3.4D, we need to discard those of 0{a''^^) in 
Eq. ( A.lD at the same time. Since the five-loop coefficient in Eq. (3.4) is not yet known we set 
to zero in our numerical analysis. 

In Fig. ^ the result for a^p{Mz) as a functions ^b is displayed for the one- to five-loop 
analysis. For illustration, fib is varied rather extremely, by almost two orders of magnitude. 
While the leading-order result exhibits a strong logarithmic behaviour, the analysis is gradually 
getting more stable as we go to higher orders. The five-loop curve is almost flat for > 1 GeV 
and demonstrates an even more stable behaviour than the four-loop analysis of Ref. [10]. It 
should be noted that around fib ^ ^ GeV both the three-, four- and five-loop curves show a 
strong variation which can be interpreted as a sign for the breakdown of perturbation theory. 
Besides the fib dependence of a^p{Mz), also its absolute normalization is significantly affected 
by the higher orders. At the central matching scale fib = Mb, we encounter a rapid convergence 
behaviour. 


4. Effective coupling between a Higgs boson and gluons 

In this Section we want to discuss the relation between and the coupling of a scalar Higgs 
boson to gluons. Due to the fact that gluons are massless, there is no coupling at tree-level. At 
one-loop order the HGG coupling is mediated via a top-quark loop depicted in Fig. 0(d). 











|a,(GeV) 


Figure 2: fib dependence of a^P{Mz) calculated from {Mr) = 0.36 and Mb = 4.7 GeV. The proce¬ 
dure is described in the text. The dotted, short-dashed, long-dashed and dash-dotted line corresponds to 
one- to four-loop running. The solid curve includes the effect of the new four-loop matching term. 


For an intermediate-mass Higgs boson which formally obeys the relation Mh 'C mt it is 
possible to construct an effective Lagrangian of the form 

Aff =- ^CiOi , (4.1) 

yyj 

with the effective operator 

Oi = {G%)^ , (4.2) 


where is the colour field strength. The coefficient function Ci incorporates the contribution 
from the top-quark loops. At one-loop order it is easy to see that the contribution from the 
triangle diagrams can be obtained through the derivative of the one-loop diagram for H^ with 
respect to the top-quark mass. However, at higher-loop orders this simple picture does not hold 
anymore and the relation between the HGG diagrams and derivatives of the two-point functions 
containing a top-quark loop gets more involved. In Ref. [10] an all-order low-energy theorem 
has been derived which establishes such a relation and which has a surprisingly simple form (for 
definiteness we specify to the top-quark in this Section): 
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An appealing feature of Eq. (^) is that at a given order in a* only the logarithmic contributions 
of C,g are needed for the calculation of Ci at the same order. Thus, from our calculation we can 
reconstruct the five-loop logarithms of Cg from lower-order terms and the j3 and 7 m functions 
governing the running of a* and the top-quark mass, respectively. This leads to the following 
result, at A^c = 3 and = \, 
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with rrit being the MS top-quark mass renormalized at the scale Note the appearance of the 
flavour-dependent part of in the five-loop contribution, whereas the corresponding coefficient 
from the anomalous mass dimension does not appear. We want to stress that the term of order 
covers the contributions from five-loop diagrams like the one in Fig. [^(e). 

Evaluating Eq. (^) numerically leads to 
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Again one observes large cancellations between the and nj term in the five-loop contribution 
to Cl. 

Note that the result of Eq. ( ^4.4|) constitutes a building block for the N^LO calculation to 
the Higgs boson production and decay in the two-gluon channel, for which the complete answer 
currently is certainly out of range. Still, the five-loop result for Ci constitutes a high-order result 
in perturbative QCD which is of theoretical interest by itself. 


5. Conclusions 


In this paper the decoupling constant of the strong coupling is presented to four-loop order. This 
constitutes a fundamental quantity of QCD and is one of the very few known to such a high 
order. The decoupling constant is necessary for performing a consistent running of ag with five- 
loop accuracy including important effects from the crossing of quark thresholds. The calculation 


has been performed analytically, and the main result can be found in Eq. (3.1). With the help 
of a low-energy theorem it is possible to derive the five-loop result for the effective coupling of 
the Higgs boson to gluons, which constitutes a building block in the corresponding production 
and decay processes. 

We want to mention that the result for in Eq. o has been obtained independently 
in Ref. [19]. Except for QGRAF, which is used for the generation of the diagrams, there is no 
common code. Even the master integrals have meanwhile been computed independently [20] 
and for the renormalization a different procedure has been chosen. 
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A. Results for 


Replacing in Eq. o the MS mass m/j by the pole mass using the three-loop approxima¬ 
tion [21-23] one gets 
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